Low-density, one-dimensional quantum gases in a split trap 
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We investigate degenerate quantum gases in one dimension trapped in a harmonic potential that 
is split in the centre by a pointlike potential. Since the single particle eigenfunctions of such a system 
are known for all strengths of the central potential, the dynamics for non-interacting fermionic gases 
and low-density, strongly interacting bosonic gases can be investigated exactly using the Fermi- 
Bose mapping theorem. We calculate the exact many-particle ground-state wave-functions for both 
particle species, investigate soliton-like solutions, and compare the bosonic system to the well-known 
physics of Bose gases described by the Gross-Pitaevskii equation. We also address the experimentally 
important questions of creation and detection of such states. 

PACS numbers: 03.75.Fi, 05.30.Jp 



I. INTRODUCTION 

Low dimensional quantum gases are fastly progressing 
from being a solely successful theoretical concept |jj LJ lil 
0, 0| towards becomingexperimentally available systems 
in the laboratories HQ, EDI • Small scale trapping 
potentials for low temperature samples allow the creation 
of degenerate fermionic or bosonic gases, for which the 
dynamics along two spatial directions can be 'frozen out'. 
Recent experimental progress has been the observation of 
quasi-one-dimensional Bose-condensates showing spatial 
phase fluctuations fll| and the generation of trapping 
frequencies in two-dimensional optical lattices of up to 
about 120kHz [IJ. 

At very low temperatures the interaction between two 
neutral particles can often be described using a hard core 
potential that is dependent on a single parameter only, 
the s-wave scattering length. The Tonks-Girardeau (TG) 
gas is a model for a one-dimensional bosonic gas made 
from particles interacting via such a hard-core potential 
d 0| and in the limit of point-like interactions it can be 
solved exactly using the Fermi-Bose mapping theorem 
H]. This theorem is one of the particularities of one- 
dimensional systems and it states that the properties of 
a strongly interacting, hard-core bosonic gas can be cal- 
culated by considering a non-interacting fermionic sys- 
tem. The first examples for which the TG gas model was 
solved were free space and box-like systems with periodic 
boundary conditions, and recently Wright and Girardeau 
have calculated its solution in an external harmonic po- 
tential [HEl 

In order to find the many-particle solutions to a given 
geometry for the TG gas (or for non-interacting fermions) 
one must know the exact single particle eigenstates. How- 
ever, since the list of exactly solvable single particle prob- 
lems in quantum mechanics is limited, there are even 
fewer many-particle problems that can be exactly solved. 
Here we present a further solvable system, which com- 
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prises of a gas trapped in a harmonic oscillator potential 
split in the center by a point-like potential. This model 
can, on the one hand, be seen as a limit of double well 
situations , which have recently been very thoroughly 
investigated for examining weakly coupled Bosc-Einstein 
condensates 0] . On the other hand it is a good approx- 
imation to the situation where an impurity is trapped in 
the center of a gas. 

In Fcrmi-Bose mapping theory the bosonic wave- 
function can be calculated directly from an appropriately 
chosen fermionic one by symmetrization. Since the sym- 
metry or antisymmetry of a wave function does not have 
an influence on the density distribution, the spatial den- 
sity profiles for fermionic and bosonic samples in this 
limit are indistinguishable. Therefore, whenever results 
concerning density distributions are presented in this pa- 
per, they apply to fermionic as well as to bosonic samples. 

The organization of the paper is the following: we will 
first briefly introduce the Fermi-Bose mapping theorem 
and its application to the TG gas. Then we will intro- 
duce the 5-split harmonic oscillator and present its exact 
solutions for arbitrary strength of the splitting potential. 
After this we will derive the exact many-particle ground- 
state wave-functions for the fermionic and the bosonic 
case in the limit of a completely split oscillator and in- 
vestigate the dynamical behaviour when the central po- 
tential strength is changed. To compare the results to 
Bose gases in the Gross-Pitaevskii limit, we investigate 
a soliton-like solution. We suggest different methods for 
experimentally realizing a system like this and outline a 
strategy to observe it. 



II. FERMI-BOSE MAPPING THEOREM 

We consider a gas of N neutral bosonic atoms of mass 
m at zero temperature in one dimension. The transverse 
dynamics is assumed to be 'frozen out' by a very stiff 
harmonic potential. The density of the gas shall be low 
enough to fulfill the requirement that the system is in the 
Tonks-Girardeau limit [TtI Il8| . In this limit the interac- 
tion between the particles can be treated as point-like. 
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Allowing that the atoms are longitudinally trapped in an 
external potential, V(x), the many-particle Hamiltonian 
can be written as 



N 
i<j 



tj\). (1) 



Here the first sum represents the single particle Hamil- 
tonian, Hq, and the second sum accounts for the inter- 
particle interaction. The coupling constant is given by 
gi = — 2h 2 /2maiu, with the repulsive, one-dimensional 
s-wave scattering length given by aiD = (— a 2 L /2a,3i))[l — 
Here a 3 D is the three-dimensional 
s-wave scattering length, a± the size of the ground 
state of the transversal harmonic trapping potential and 
C = 1.4603 .... In the Tonks-Girardeau limit the one- 
dimensional coupling between the particles becomes large 
and we will consider g\ — ► oo in the Hamiltonian l(T)l. In 
this situation, instead of including the scattering term 
into the Hamiltonian, it is more convenient to account 
for the strong repulsive interactions using a constraint 
on the allowed solutions 1141 



ip = 



if 



= o, 



(2) 



A sketch of the potential is shown in the upper left of 
Fig. ^ One can note immediately that the odd eigen- 
functions, j 02n+i, of the undisturbed harmonic oscilla- 
tor (k = 0) are still eigenfunctions of the Hamilto- 
nian Ho, since they are not affected by the potential 
(V-'2ri+i(0) = 0). In particular, they are independent of 
the strength k of the 5-function. 

To find the even eigenfunctions of the Hamiltonian in 
eq. Q, let us scale all quantities in units of the ground 
state size ao = ^/h/2muj and the energy eo = Tko of 
the harmonic oscillator with k = 0. The stationary 
Schrodinger equation then takes the form 



dx 2 



1 



k5(x) 



<f>n(x) = 0, 



(5) 



where we have also introduced the normalized energy 
e n = —E n /(%uj) and the normalized delta function 
strength k — Kao/eo- Note the negative sign of the nor- 
malized energy, which has been chosen for easier mathe- 
matical treatment. One can see that for x > eq. (J5J is 
equal to the differential equation for the parabolic cylin- 
der functions and the solutions that have the physically 
correct asymptotic behavior are Whittaker's functions, 
U(e n ,x), given by [2Cj 



Since this constraint is equivalent to the demand of the 
Pauli exclusion principle for a gas of spin- less fermions, 
one can calculate the bosonic eigenfunctions, ipB, of the 
Hamiltonian of eq. £Q) from the fermionic eigenfunctions, 
ipFt of the- non- interacting Hamiltonian by proper sym- 
metrisation 



*b(xi, . . .,x N ) = |*_f(xi, . . . ,Xjv)|. 



(3) 



This procedure is known as the Fermi-Bose mapping the- 
orem 0, 0|. It can be generalized to time dependent 
problems. 

To find the many-particle wave-function of the inter- 
acting Bose-gas we therefore need to solve the Slater de- 
terminant for non-interacting fermions. This determi- 
nant is composed of the eigenstates of the corresponding 
single particle Hamiltonian, Hq- However, the number of 
physical systems where the single particle eigenstates are 
available analytically is very small. One example is the 
(5-split harmonic oscillator and we will review its single 
particle states in the next section. 
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The M(a, 6, c) are the confluent hypergeometric func- 
tions. To first construct the symmetric eigenfunctions 
of the Hamiltonian in eq. on the domain x ^ we 
can now define [jjl 



(f> n (x) = CU(e n , \x\), 



(9) 



with the normalization constant C chosen such that the 
wave-function is normalized to one. Finally, evaluating 
the continuity condition at <j) n {x = 0) 



d 

dx 



0n(O+) 



d 

dx 



<j> n (0-) = K0„(O), 



(10) 



we get an implicit relation between the strength of the 
((-function and the eigenvalues of the Hamiltonian 



III. EIGENSTATES OF THE 5— SPLIT 
HARMONIC TRAP 

The single particle Hamiltonian describing a one- 
dimensional harmonic trap of frequency lo with a 8- 
potential of strength k in the center is given by 



Ho = -^ 2 +lm^ + K 6(x). (4) 



r f2 + i e 
r(| + ie 



(11) 



Since the T-function has poles for negative integer values, 
one can immediately sec that the energy eigenvalues for 
a strong central barrier (n — > oo) converge to E n — > 

§?iw, Qfou), They become identical to the values 

of the odd eigenfunctions in the limit At = oo, leading 
to a double degeneracy of all eigenstates 0. Fig. Q] 



3 



K 



P0.3 




p0.3 




linear operations, the determinant can be reduced to 
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FIG. 1: Sketch of the delta split harmonic oscillator and the 
respective single particle ground state (n = 0) densities for 
different values of k. The functions for n — 0, 1, oo corre- 
spond to the energy eigenvalues of E n — \TijU), 1.087jw, §7i6l> 
respectively. 



displays the square of the ground state eigenfunctions of 
the Hamiltonian in eq. I@J for different values of k. One 
can clearly see the development of a node at x = for 
increasing strength of the splitting potential. 

In the limit k = oo, when the neighbouring even and 
odd eigenfunctions become degenerate, the Whittaker 
functions of the even states simplify to Hermite poly- 
nomials with an absolute value argument. 



IV. QUANTUM GAS IN A 5-SPLIT 
HARMONICS TRAP 



To find the exact many-particle wave-function for 
bosons or fermions, one has to calculate the Slater de- 
terminant from the single particle eigenfunctions 4> n {x) 
and, for the bosons, apply the symmetrization procedure 
to it. Due to the complicated structure of the Whittaker 
functions for finite values of n this can only be done nu- 
merically. However, for the limiting case of an infinitively 
high central (^-function this can in fact be done analyti- 
cally. The eigenfunctions are in this case given by 



x 2 

ip n (x) = C n e~~ H n (x) 
ipn(x) = C n+1 e 2 H n+1 (\x\) 



for n odd, (12a) 
for n even, (12b) 



with the normalization C n — (y / 7rao2"n!) 2 ) and the 
Slater determinant is defined by 



tpF{xi, . . . ,xn) 



1 



(N-l,N) 

det ijj n (xj). 



'N\ (nj) = (0,l) 



(13) 



Inserting the explicit expressions for the Hermite poly- 
nomials into the determinant and applying basic multi- 



(N-i,N) t 
det ib n (xj) — C 2 
(n,j)=(0,l) 



\Xl\ Xi \Xi\ 



(14) 



where C = 11^=1 Cfn-i" Unfortunately in this form the 
determinant cannot be easily solved, due to the appear- 
ance of the absolute values. 

We can, however, rewrite the problem in a more 
amenable form. It is easy to see that for the case of 
k = 00 due to the double degeneracy of each energy 
level, one can choose a new set of eigenfunctions given by 
(f>n = Tffen ±i>2n+i) .15] ■ These eigenfunction corre- 
spond to the particles being physically located entirely ei- 
ther left or right of the barrier and the wave- functions (jx^ 
are mutually symmetric with respect to x — 0. Therefore 
we can restrict our treatment to the spatial area x > 
and consider only ^ particles (for simplicity we assume 
an even number of particles). The Slater determinant 
can then be explicitly calculated by just employing the 
odd Hermite polynomials of eq. (|12afl 



(f-i,f) 

det H 2n +i{xj) 
(«j)=(o,i) 
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which immediately gives a simple and exact expres- 
sion for the fermionic many-particle wave- functions of 
eq. (I13fl . The bosonic form then follows directly from 
eq. (0 and is given by 



* B (xi, . . . ,x N ) 
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Co 



U,k) = (l,j+1) 



It resembles the wave-function 



with C = n„=i ^2n-l. 
found by Girardeau et al. [2 for the undisturbed har- 
monic oscillator, and differs due to the reason that we 
have constructed the Slater determinant only from the 
odd Hermite polynomials. 

In Fig. |2 the density of a cloud of N = 50 bosons is 
shown for different heights of the central potential. One 
can clearly see that the central notch gets deeper with 
increasing values of k, but one also finds an increase in 
density of about 25% just around the hole for large values 
of k. As one would expect, the influence of the (5-potential 
only extends a short distance (the 'healing length') from 
x = and does not have any influence on, say, the size 
of the atomic cloud. 
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FIG. 2: Ground state density for a one dimensional quantum 
gas of 50 particles in a 5-split harmonic trap. The strength of 
the 5-function is given by re. Note the increase in density in 
the center of the trap around the notch. Note also that the 
appearance of the 5-function has no influence on the size of 
the atomic cloud. 



The density at x = of the g function of re is 

displayed in the LHS of Fig.[3J One can see that for values 
of re of a few to the density drops quickly towards zero. 
On the RHS the full-width-half-max size of the central 
notch is shown. As compared to the solitonic solutions 
for Bose gases described by the Gross-Pitaevskii equation 
where the width of the notch decreases with increasing 
darkness, the opposite behaviour is found here. This is 
of course due to the fact that the density void of the 
soliton in a Bose condensate is due to a phase jump in 
the wave-function and not an external potential. 



V. DARK SOLITONS 

As already mentioned in the last section, the density 
distributions shown in Fig.[21for finite re resemble the dark 
and grey soliton states of the Gross-Pitaevskii equation 
(GPE). Dark (and grey) solitons are solutions to the GPE 
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FIG. 3: The left hand side shows the density of a gas at x — 
depending on the strength of the delta function. The right 
hand side shows the FWHM of the central notch. Note that 
the width increase with increasing 'darkness'. 



with repulsive interaction and their existence in dilute 
Bose-condensates has been experimentally demonstrated 
pll El ] . Recently also the appearance of dark solitons in 
TG gas has been examined in a toroidal trap |2^| . 

Having derived the exact many-particle state for the 
pierced harmonic oscillator, we can now compare the 
TG 'notch' state with soliton solutions of the Gross- 
Pitaevskii equation in a trap. Note, that we have not 
only the 'dark notch' solution at hand, but we also know 
the full family of solutions for 'grey notches', where the 
density is reduced but does not go to zero. 

We are, however, again not able to derive any analyt- 
ical expression for the grey soliton densities, since the 
summations over the eigenstate densities involve a sum- 
mation over e„. This doesn't represent a series of integers 
nor is there an easy relationship between the integer n 
and the e n and we therefore resort to numerical simula- 
tions. 

To simulate the time evolution of the bosonic gas, one 
can make use of a time dependent version of the Fermi- 
Bose mapping theorem . The many-particle fermionic 
wave-function is then given by 



tp F (xi, ■ ■ - ,x N ;t) 



(N-l.N) 

C det <f>i{xj,t), (17) 

(i,j)=(0,l) 



and the single particle eigenstates of the external poten- 
tial evolve under the time-dependent Schrodinger equa- 
tion 



ih 



d(j)j{x, t) 

at 



IL^L v(x t) 

2m dx 2 ' 



<t>i{x,t). (18) 



Symmetrization of these states then leads to the appro- 
priate solutions for the bosonic case. To investigate the 
dynamics of the notch, we assume that we can prepare 
the system in the stationary state we want for fixed re 
and then we instantaneously switch re = 0. 

First, it is easy to understand that the recurrence fre- 
quency for the density notch in the gas is given by the 
trap frequency. Since the density distribution is in partic- 
ular the density distribution for non-interacting fermions, 
after one trapping period lo every fermion is back in the 
same state where it started and the notch is recreated 
(see Fig. QJ. Even though this result is to be expected 
for non- interacting fermions, it is not at all obvious for a 
system of strongly interacting bosons. In particular it is 
in contrast to behaviour of dark solitons in the GPE, for 
which the oscillation frequency is given by lo/\/2 j2J]. 

Since the filling in and the oscillation of the dark notch 
as seen in Fig.0]are purely due to the even eigenfunctions 
(the odd ones are stationary states for any value of re), 
one could consider to stabilize the density notch by im- 
printing a phase-step with A<f> = 7r onto the sample. This 
would artificially create the same stabilization mecha- 
nism that the soliton solution of the GPE has, i.e. it 
would transform the single particle states with an even 
parity to states with an odd one. However, since one 
cannot make this process parity selective, at the same 
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FIG. 4: Density at x = during one oscillation period after 
re has instantaneously been changed from re = oo to k = 0. 
One can clearly see the recurrence of the notch after one trap 
period. 



time it would transfer the odd states into even states 
and nothing would be won. 



VI. DETECTION 

The absence of a self-similar free expansion has been 
pointed out as one criterion to distinguish the TG gas 
from a one-dimensional Bose-gas in the Gross-Pitaevskii 
regime [2£j. Since the atomic cloud is only disturbed on 
a very small length scale by the central 5-function, and, 
as can be seen from Fig. the lifetime for the density 
notch once the central barrier is removed is much shorter 
than a trap frequency, the modification due to the split 
trap has no measurable effect on the size and dynamics 
of the expanding cloud (see Fig. |5j. Especially, since 
the lifetime of the density notch is so short, once cannot 
observe any broadening of it 

This fact can also be seen by realizing that the energy 
increase for finite re is at maximum AE — > ^Tiuj for the 
limiting case of k — > oo, which is very small on the scale 
of the total energy of the cloud, E to t = N 2 /2. (Note that 
the chemical potential alone is already [i — (N —l/2)Tno.) 

We suggest here a different approach that relies on the 
insensibility of the state towards phase imprinting. In 
standard traps (re = 0) phase imprinting has a substan- 
tial effect on the density distribution of a sample, a fact 
which has been successfully utilized for creating dark soli- 
tons in ground state Bose-Einstein condensates. In the 
dynamics following the imprinting the gas tries to ad- 
just its density distribution to the phase, thereby excit- 
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FIG. 5: Integrated absolute difference between the density 
distributions of the ground state gas for re = and re = oo 
during the free expansion. The fast decay is a sign of the 
indistinguishability of the two ground states via the density 
distributions. The inset shows the density p(x — 0) for re = oo 
(broken line) and for re = (full line). One can see that the 
hole vanishes at very early stages during the free evolution. 
Due to the absence of an external potential, there are no 
recurrences here. The timescale in both graphs is in units of 
the inverse trap frequency. 



ing various phonon modes and modifying the Thomas- 
Fermi density distribution. Recently phase imprinting 
was theoretically investigated for fermionic samples, and 
the emergence of solitonic structures was also found there 
|26j. However, since the phase evolution of the single par- 
ticles in a Fermi sea is not coherent, the effects found were 
less pronounced. 

To demonstrate the insensibility of the bosonic and 
fermionic states in a <5-split trap, we consider here the 
limiting case of re — > oo, and the reaction to an imprinted 
phase of Acf> = tt. Applying this imprinting only spatially 
selectively for, say, x < 0, one can see immediately that 
the eigenfunctions ip n {x) with odd values of n turn into 
VVi-i(M) and vice versa. Since after the imprinting we 
again deal only with stationary eigenfunctions, the phase 
imprinting (assuming an ideal imprinting mask) does not 
have any observable effect on the density distribution. 
Mathematically one can see immediately that the Slater 
determinant due to its multilinearity is unchanged by the 
imprinting of the phase (for an even number of particles) . 
Therefore the absence of distortions could be used a a 
possible tool to detect the split trap state. 



VII. EXPERIMENTAL REALIZATION OF THE 
(5-POTENTIAL 

Let us finally remark on the experimental realization of 
such a potential, since the <5-function is usually regarded 
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as a mathematical idealization far from experimental re- 
ality. A very successful method to manipulate the den- 
sity profile of an atomic cloud is to employ dipole forces 
with appropriately detuned lasers. It is nevertheless ex- 
perimentally quite challenging to focus and stabilize a 
laser beam to such a small diameter that it represents 
a point-like potential to the atoms as necessary for the 
situation considered here. However, since low densities 
are best reached in shallow traps, the ratio of the size 
of the single particle states vs. the fwhm of the focused 
laser beam increases with increasing trap size, thereby 
creating a more and more localized potential. The tun- 
ing of k could in this case be done 'easily' by adjusting 
the intensity of the laser beam. 

Here we suggest an alternative method which at the 
same time also allows for easy tuning of k, even though 
it is also experimentally more challenging. 

If optical potentials do not suffice in terms of localiza- 
tion, the obvious next step is to employ potentials made 
from matter waves. Multicomponent condensates have 
been investigated in the last years and the conditions for 
such systems to exhibit phase separation have been de- 
termined. Depending on the internal energy of the atoms 
and on their respective masses shell like structures have 
been predicted, where one component occupies the inner 
part of the trap and the other component occupies the 
outer part. By realizing a situation like this and then re- 
ducing the number of atoms occupying the central part of 
the trap, one can, in the ideal case, come to the situation 
where a single atom is trapped in the trap center and 
the second component is arranged around it. This would 
correspond to an almost perfect point-like potential for 
the atoms of the outer component. With the help of a 
selective Feshbach resonance the strength of the interac- 
tion between the single atom and the gas could then be 
adjusted. 

VIII. CONCLUSIONS 

We have investigated the behaviour of low density, 
one-dimensional Bose and Fermi gases in a <5-split har- 



monic trap. We have derived the exact many-particle 
wave function for non-interacting fermions and were able, 
due to the Fermi-Bose mapping theorem, to determine 
also the exact expression for the many-particle wave- 
function for the bosonic TG gas in the limiting case of 
an infinitely high <5-barrier. Both many-particle wave- 
functions turned out to be simple products over quadratic 
terms. We have thereby presented another of the rare ex- 
amples of exactly solvable models for many-particle sys- 
tems. 

We have further numerically investigated the dynami- 
cal behavior of many-particle, split-trap eigenstate when 
the central barrier is suddenly removed and compared it 
to the dynamics of the dark soliton solutions in the GPE. 
The recurrence time for the 'notch' in the linear low den- 
sity limit was found to be the trapping frequency, which 
is different to the oscillation frequency of the dark Gross- 
Pitaevskii soliton of uj/^/2. 

Finally, we have suggested using the systems insen- 
sitivity to phase imprinting as a tool for detecting the 
presence of the notch-state and employing a second com- 
ponent and a selective Feshbach resonance to create and 
manipulate the point-like central potential. 
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